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(54) Title: METHODS AND APPARATUS FOR BLIND SIGNAL SEPARATION 

(57) Abstract 

A set of generalized architectures, frameworks, algo- 
rithms, and devices for separating, discriminating, and recov- 
ering original signal sources by processing a set of received 
mixtures and functions of said signals based on processing of 
die received, measured, recorded or otherwise stored signals 
or functions thereof. There are multiple criteria that can be 
used alone or in conjunction with other criteria for achiev- 
ing the separation and recovery of the original signal content 
from the signal mixtures* The system of the invention en- 
ables the adaptive blind separation and recovery of several 
unknown signals mixed together in changing interference en- 
vironments with very minimal assumption on the original sig- 
nals. The system of this invention has practical applications to 
non-multiplexed media sharing, adaptive interferer rejection, 
acoustic sensors, acoustic diagnostics, medical diagnostics and 
instrumentation, speech, voice, language recognition and pro- 
cessing, wired and wireless modulated communication signal 
receivers, and cellular communications. 
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METHOPS AN P APPARATUS FOR 
PLIWP SIGNAL SEPARATION 

CROSS-REFERENCE TO RELATED APPLICATION 
This application is a continuation-in-part of the U.S. provisional patent 
application entitled ARCHITECTURES, FRAMEWORKS, ALGORITHMS AND 
DEVICES FOR REAL TIME ADAPTIVE SIGNAL SEPARATION, 
DISCRIMINATION AND RECOVERY, Serial No. 60/050,147, filed June 18, 1997, 
invented by Gamze Erten and Fathi M. Salam. 

BACKGROUND OF THE INVENTION 

Field of the favertfjon 

The present invention relates to methods and apparatus for recovering original 
signal information or content by processing multiple measurements of a set of mixed 
signals. More specifically, the invention relates to adaptive systems for recovering 
several original signals from received measurements of their mixtures. 

Discussion of Related Art 

The recovery and separation of independent sources is a classic but difficult 
signal processing problem. The problem is complicated by the fact that in many 
practical situations, many relevant characteristics of both the signal sources and the 
- mixing media are unknown. 

To best understand the problem, the following problem statement is helpful: 
With reference to FIGURE 1 of the drawings, consider N independent signals 102, 
s,(t), ... , and s N (t). Independent signals 102 may represent any of, or a combination 
of, independent speakers or speeches, sounds, music, radio-based or light based 
wireless transmissions, electronic or optic communication signals, still images, 
videos, etc. Independent signals 102 may be delayed and superimposed with one 
another by means of natural or synthetic mixing in the medium or environment 104 
through which they propagate. A signal separation process 106 transforms the 
mixture of independent signals 102 into output signals 108, u,(t), ... , and u N (t). 
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Two main categories of methods are often used for the recovery and separation 
of independent sources: neurally inspired adaptive algorithms and conventional 
discrete signal processing. 

Neurally inspired adaptive architectures and algorithms follow a method 
5 originally proposed by J. Herault and C. Jutten, now called the Herault-Jutten (or HJ) 
algorithm. The suitability of this set of methods for CMOS integration have been 
recognized. However, the standard HJ algorithm is at best heuristic with suggested 
adaptation laws that have been shown to work mainly in special circumstances. The 
theory and analysis of prior work pertaining to the HJ algorithm are still not sufficient 

10 to support or guarantee the success encountered in experimental simulations. Herault 
and Jutten recognize these analytical deficiencies and they describe additional 
problems to be solved. Their proposed algorithm assumes a linear medium and 
filtering or no delays. Specifically, the original signals are assumed to be transferred 
by the medium via a matrix of unknown but constant coefficients. To summarize, the 

15 Herault-Jutten method (i) is restricted to the full rank and linear static mixing 

environments, (ii) requires matrix inversion operations, and (iii) does not take into 
account the presence of signal delays. In many practical applications, however, 
filtering and relative delays do occur. Accordingly, these methods fail to successfully 
separate signals in many practical situations and real world applications. 

20 Conventional signal processing approaches to signal separation originate 

mostly in the discrete domain in the spirit of traditional digital signal processing 
methods and use the statistical properties of signals. Such signal separation methods 
employ computations that involve mostly discrete signal transforms and 
filter/transform function inversion. Statistical properties of the signals in the form of a 

25 set of cumulants are used to achieve separation of mixed signals where these 

cumulants are mathematically forced to approach zero. This constitutes the crux of 
the family of algorithms that search for the parameters of transfer functions that 
recover and separate the signals from one another. Calculating all possible cumulants, 
on the other hand, would be impractical and too time consuming for real time 

30 implementation. 

The specifics of these two methods are elaborated on below. 
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NemaHy inspired archijeyfrire frfld algorithms for signal gepflratipn 
These sets of neurally inspired adaptive approaches to signal separation 
assume that the "statistically independent" signal 

vector £(0 = [j,(0» ... , and s N (t)] T is mixed to produce the signal vector M(t) . 
5 The vector M(t) is received by the sensors (e.g. microphones, antenna, etc.). 

Let the mixing environment be represented by the general (static or dynamic) 
operator!. Then, 

M(t) = t(S(0) Equation (1) 

There are several formulations that can be used to invert the mixing process, i.e. 
10 operator x, in a "blind" fashion where no a priori knowledge exists as to the nature or 

content of the mixing operator t or the original sources S(t) . We group these into 

two categories, static and dynamic. Additional distinctions can be made as to the 

nature of the employed adaptation criteria, e.g. information maximization, 

minimization of high order cumulants, etc. 
15 The static case. The static case is limited to mixing by a constant nonsingular 

matrix. Let us assume that the "statistically independent" signal vector 
£(0 = [*i(0» •» > and s N i t )] T is mixed to produce the signal vector M(t) . 

Specifically, let the mixing operator x be represented by a constant matrix A , 

namely 

20 Af(r) = AS(t) Equation (2) 

FIGURES 2A and 2B show two architectures of the signal separation and 
recovery network in case of static mixing by a mixing matrix A . U(t) is the output 
which approximates the original source signals S(t) . Y(t) contains the values that 
are used in updating the parameters of the unmixing processes, i.e. W in FIGURE 

25 2Aand D in FIGURE 2B. The architecture in FIGURE 2A necessarily computes 
the inverse of the constant mixing matrix A , which requires that A is invertible, 
i.e. A ~* exists. The architecture in FIGURE 2B does not impose the restriction in 
that upon convergence the off diagonal elements of the matrix D are exactly those 
of the off diagonal elements of the matrix A . In this case, however, diagonal 

30 elements of the matrix A are restricted to equal " 1 .0." By setting the diagonal 
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elements of D to zero, one essentially concludes that the mixing process is 
invertible even if the mixing matrix is not. As shown in FIGURES 2A and 2B, the 
weight update utilizes a function of the output U(t) . The update of the entries of 
these two matrices is defined by the criteria used for signal separation, discrimination 
5 or recovery, e.g. information maximization, minimization of higher order cumulants, 
etc. 

As an example, one possible weight update rule for the case where 
U(t) = WM(t) Equation (3) 

could be 

10 w.. r = t\[W~ t + Z^M\ Equation (4) 

*'(«) 

where r\ is sufficiently small, g is an odd function, and M is the set of 
mixtures, U is the set of outputs which estimate the source signals, subscript T 
denotes transpose, and -T denotes inverse of transpose. Note that the 
function g plays an additional role in the update which can be related to the above 
15 diagram as 

HCO = g(U(t)) Equation (5) 

One uses Equation (4) to update the entries of W in Equation (3). Through this is 
an iterative update procedure, the entries of W converge so that the product WA is 
nearly equal to the identity matrix or a permutation of the identity matrix. 
20 On the other hand, in FIGURE 2B, one potentially useful rule for the update of 

the D matrix entries is generically described as 

<t (J = n/tyto)^/')) Equation (6) 

where n is sufficiently small. In practice some useful functions for / include a cubic 
function, and for g include a hyperbolic tangent function. When using this 
25 procedure, one computationally solves for U(t) from Equation (7) below 

U(t) = [/ + D]' l M(t) Equation (7) 
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at each successive step and sample point. This computation is a potentially heavy 
burden, especially for high dimensional D . 

The dynamic case. The dynamic mixing model accounts for more realistic 
mixing environments, defines such environment models and develops an update law 
to recover the original signals within this framework. 

In the dynamic case, the matrix A is no longer a constant matrix. In reference 
to the feedback structure of the static example, it is simpler to view Equation (7) 
where U(t) = [J + D]' l M(t) as an equation of the fast dynamic equation 

zU(t) = -U(t) - DU(t) + M(t) Equation (8) 

This facilitates the computation by initializing the differential equation in Equation (8) 
from an arbitrary guess. It is important however to ensure the separation of time 
scales between Equations (8) and the update procedure like the one defined by 
Equation (6). This may be ensured by making r\ in Equation (6) and x in Equation (8) 
sufficiently small. 

If we assume the dimensionality of M(t) is N, a set of differential equations 
that define the dynamic signal separation algorithm can be written as 



j+t Equation (9) 

for i = 1,...,AT 



This enumerates N differential equations. In addition, the adaptation process for the 
entries of the matrix D can be defined by multiple criteria, e.g. the evaluation of 
functions / and g in Equation (6). FIGURE 3 is a pictorial illustration of the 
dynamic model in feedback configuration. £/(/) approximates S(t) . The 
Sanction g defines the criteria used for weight update of the feedback network. 

Application of signal separation criteria for adaptation of parameters. Little 
has been outlined in the way of procedures for the application of adaptation criteria 
within the architectures defined thus far. Two implied procedures have been noted: 

First is the application of the signal separation functions, adaptation 
procedures and criteria to arbitrary points of data - regardless of whether each of these 
points is practically and physically accessible or not. Thus, the adaptive separation 
procedure applies the adaptation functions and criteria to each element of the 
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measured mixed signals individually and instantaneously, after which appropriate 
parameter updates are made. 

The second type of procedure described in FIGURE 2 A uses Equation (3). In 
this case, the criteria is applied to the entire data set, or selected data points from the 
5 entire data set. Thus, the related adaptation process does not progress per sample, but 
utilizes the whole data set over which a constant, static mixing matrix is assumed to 
apply. Although this method is somewhat more robust than the first, it is essentially 
an offline method not suitable for real time signal separation. Furthermore, when the 
assumption of a static constant matrix is incorrect, the accuracy of the unmixing 
10 process suffers. 

The transfer function based approach to signal separation 

The representation of signal mixing and separation by transfer functions makes 
this approach a dynamic environment model and method. 

A structure for separating two signals by processing two mixture 
15 measurements is illustrated in FIGURES 4A and 4B. FIGURES 4A and 4B show a 

conventional transfer function representation for a signal mixing and separation for a 
two signal system. The two signals U l and U 2 approximate S x and S 2 . G 
inverts the mixing process modeled as H . Such a transfer function approach is 
neither practical or extendible in the case of higher dimensional signals greater than 
20 two dimensions. Furthermore, the extension of the mixing environment to the transfer 
function domain has also eliminated the time domain nature of the signals. This also 
causes the exclusion of the initial conditions from the set of equations. 

These architectures for the separation functions in the transfer function domain 
result in three serious shortfalls which are all impediments to the design and 
25 implementation of a practical method and apparatus. First, this formulation, as 
expressed, precludes the generalization of the separation procedure to higher 
dimensions, where the dimensionality of the problem exceeds two. In other words, a 
practical formalization of the separation method does not exist when there are more 
than two mixtures and two sources. One can illustrate this by direct reference to 
30 current methods, where matrix multiplication terms are written out, so that each scalar 
equation defines one of the entries of the resulting product matrix desired to be equal 
to zero. Since permutations of a diagonal matrix are also allowed, multiple sets of 
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equations are created. For a two mixture problem, this results in two pairs (four total) 
of equations, each with two product terms. Beyond that the number of equations 
increases. To be precise the number of equations needed to describe the number of 
equations for a specific permutation of the N dimensional case is equal to (N 2 - N). 
5 For the two dimensional problem this value is 2. 

Second, the inversion procedure for the transfer function is ad hoc and no 
recipe or teaching exists. The impact of dimensionality plays a crucial role in this. It 
is apparent from the method that the resulting architecture gives rise to networks 
requiring transfer components whose order is dependent on products of the transfer 

10 components of the mixing environment. Thus, one cannot design a network 
architecture with a fixed order. 

Third, the initial conditions cannot be defined since the formulation is not in 
the time domain and cannot be initialized with arbitrary initial conditions. Hence, the 
method is not suitable for real time or on line signal separation. 

15 What is needed is a method to separate mixed signals through media or 

channel wherein a high quality of signal separation is achieved. What is needed is a 
method to recover and separate mixed signals transmitted through various media 
wherein the separation of signals is of such high quality as to substantially increase 
(i) the signal carrying capacity of the medium or channel, (ii) the quality of the 

20 received signal, or (iii) both. The media or channels may include a combination of 
wires, cables, fiber optics, wireless radio or light based frequencies or bands, as well 
as a combination of solid, liquid, gas particles, or vacuum. 

SUMMARY OF THE INVENTION 
25 The present invention is directed toward a signal processing system for 

separating a plurality of input signals into a plurality of output signals. The input 
signals are composed of a mixture of a plurality of source signals that have been 
affected by a medium. The source signals are associated with a plurality of sources. 
The output signals estimate the source signals. The system comprises a plurality of 
30 sensors for detecting the input signals and a storage device for receiving and storing 
the input signals. The system also comprises an architecture processor for defining 
and computing a signal separation architecture, where the signal separation 
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architecture defines a relationship between the input signals and the output signals, 
and the relationship has both constant parameters and time varying parameters. The 
system also comprises an update processor for computing a rate of change for each 
time varying parameter and the time varying parameters in response to the rate of 
change associated with each time varying parameter, and an output processor for 
computing the output signals based on the signal separation architecture, the constant 
parameters, and the time varying parameters. 

BR IEF PESCWnON OF TH E DRAWINGS 
FIGURE 1 shows a graphical depiction of the signal separation, discrimination 

and recovery problem. 

FIGURES 2A and 2B show two static neural network structures for signal 

separation. 

FIGURE 3 shows an architecture of the signal separation and recovery 
network in case of feedback dynamic mixing and separation models. 

FIGURES 4A and 4B show a conventional transfer function representation 
and method for signal mixing and separation for a two signal system. 

FIGURE 5 A shows a process flow diagram of a method of the present 
invention. 

FIGURE 5B shows a conventional transfer function frequency domain 
representation for signal mixing and separation. 

FIGURES 6A and 6B show two mixing models for the state space time 
domain architecture. 

FIGURES 7A and 7B show two signal separation models for the state space 
time domain architecture. 

FIGURE 8 shows a typical DSP implementation architecture. 

FIGURE 9 shows a typical DSP internal architecture. 

FIGURE 10 shows multithreaded adaptation processes. 

FIGURE 1 1 shows one embodiment of a constant width long stack based on a 
first in first out (FIFO) stack construct and data acquisition process. 

FIGURE 12 shows another embodiment of a constant width long stack based 
on a first in first out (FIFO) stack construct. 



8 



WO 98/58450 



PCT/US98/13051 



FIGURES 13A and 13B show an application of the separation criteria using a 
set of data points (short stack) from a long stack which contains all relevant sets of 
data points, measured or computed over a window of time. 

FIGURE 14 shows an audio application based on the signal separation and 
5 recovery procedures of this invention. 

FIGURE 15 shows three audio input interfaces to the DSP based handheld 
device that can be embedded into the microphone-DSP interface. 

FIGURE 16 shows a schematic of an audio device suitable for use with the 
present invention. 
10 FIGURE 1 7 shows a smart microphone array device. 

DETAILED DESCRIPTION OF THE INVENTION 
FIGURE 5 A shows a process flow diagram of a method of the present 
invention. This includes: (1) obtaining samples (block 502); (2) interpolating or 

15 upsampling (block 504); (3) selecting sample set (block 506); (4) computing adaptive 
parameters (block 508); (5) computing internal states (block 510); (6) computing 
outputs (block 512); and (7) storing and/or presenting outputs (block 514). 

Obtaining samples includes obtaining the multi-channel data recorded through 
multiple sensors, e.g. microphones. Such data could also come from previously 

20 mixed sound tracks, e.g. recording studio mixed signals. Data can be sampled on line 
for a real time process or recalled from a storage medium, tape, hard disk drive, etc. 

Interpolating is performed for algorithms that may require oversampling. If 
oversampling of the signal is not possible (e.g. prerecorded sound tracks) or not 
practical, one can interpolate between sample points to come up with more data points 

25 for the same duration of the signal. Upsampling can be performed depending on the 
sample rate of the obtained data and the target input sample rate, by upsampling by 
certain factor. For example, for making up data of 60 kHz from data sampled at 10 
kHz, the upsampling factor is 6. In other words, five points are generated in between 
each sample point available. 

30 Selecting sample set involves selecting the data points used at each iteration of 

the process. 
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Computing adaptive parameters may involve a method which uses the 
derivatives of a function to compute the value of the function (the function is the 
adaptive parameters in this case). Most of these methods are referred to as integration 
methods for solving differential equations, e.g. Runge Kutta. 

Computing internal states involves looking at the architecture. These internal 
states may also be in the form of actual states and derivatives of these states, or in the 
form of samples in time. State values can be computed from their derivatives using 
various types of integration methods for solving differential equations, e.g. Runge 
Kutta. 

Computing outputs uses the states and parameters computed earlier. 

A distinction can be made in categorizing the architectures of the present 
invention. This pertains to the mixing and the separation process models and defines 
the procedure of separation discrimination or recovery - depending on the application 
environment and nature of the measured signals. This distinction sets the 
mathematical equations and axioms about both the mixing and separation models and 
thus applies to both the mixing and the separation procedures. 

There are three sets of architectures and methods which can be itemized as 
follows: 

1 . Transfer function-based frequency domain architecture and method. 

2. State space time domain architecture and method: 

3. Mapping architecture and method. 

Transfer function based frequency domain signal separation architecture and 
method 

FIGURE 5B shows a conventional transfer function frequency domain 
representation for signal mixing and separation. U(s) and S(s) are 
multidimensional signals, and U(s) approximates S(s) . H(s) inverts the mixing 
process which is modeled as H(s) . The computation of this inverse may involve 
construction of conventional signal filters or maximization or minimization functions 
or criteria based upon optimal control theory and calculus of variations. 

For the relationship between the measured mixed signals M(t) and the 
original sources of signals S(t) to be represented as a transfer function, one must 
assume that the environment is linear time invariant. This is because one can 
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represent only linear time invariant relationships by a transfer function. Thus, 
formulation of the signal separation problem by a transfer function, while including 
dynamics and filtering effects, also assumes time invariant (stationary) environment. 

There are, however, many applications of this formulation to well defined, 
band-limited signals about which many such assumptions can be made. 

The relationship between the Laplace transform of the source signals S(t) , 
i.e. S(s) and Laplace transform of the measured mixed signals M(t) , i.e. M(s) 
are defined as follows: 

= H(s) Equation (10) 

S(s) 

The objective of this invention is, without knowing the actual H(s) , to 
obtain H(s) as the inverse of H(s) and thus reconstruct the original signals. 

The execution of this process for two signals, i.e. where the dimensionality of 
both S and M is two, is possible because for the product 

1 0 

H(s)H(s) = ^ Equation (11) 



or 



H(s)H(s) = 



0 1 

1 0 



Equation (12) 



For this two dimensional case, the set of equations that need to be satisfied are 
easily solvable. In extending the architecture to more than two signals one equates: 

H(s)H(s) = T(s) Equation (13) 

where T(s) is of diagonal form. If one equates T(s) to the identity, where all 
diagonal entries equal 1 .0, then 

H(s) = H(s) x Equation (14) 

and therefore 
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H(s) = \z adj(H(s)) Equation (15) 

det(i/(s)) 



where adj(H(s)) is the matrix with (ij)th entry computed as the (j>i)th cofactor. In 
turn, the 0,i)th cofactor is (-iy +i the determinant of the matrix formed by removing the 
jth row and the ith column of H(s) . One can further define 

Equation (16) 



Thus, a generalized structure for higher dimension emerges as a network 
defined by Equation (16). Note that, in this case, one has: 

H(s)H(s) = det (#(*))/ Equation (17) 

The actual formulation of the elements h^s) of the transfer function H(s) is 
dependent upon the characteristics of the mixing medium and the behavior and 
interactions of the original signal sources or their components. 

The processes for the construction of the elements of the separation 
filter H(s) would require that the product H(s)H(s) be an identity matrix or a 
permutation thereof, or alternately a diagonal matrix as shown in Equation (1 7). The 
computation of this inverse becomes cumbersome beyond the two dimensional case. 
The addition of each signal creates a substantial computational burden. Nevertheless, 
if one knows the specifics of the signals and the mixing medium this procedure may in 
theory be followed. 

An example procedure for the utilization of this architecture for signal 
separation discrimination and recovery can be summarized as follows: 

1 . Construct and allocate a memory structure which is used for storing the 
measurements as well as the results of the data obtained from 
processing of the measurements. 

2. Construct a transfer function N x N element matrix to be 

the H(s) transfer function (inverse of the mixing transfer function), 
the entries of which are the scalar transfer functions. 
3 Receive a new set of measurements, or mixed signals. 
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4. Process the new set of measurements. 

5. Evaluate signal separation criteria and functions. 

6. Update filter and transfer function element parameters accordingly. 

7. Accumulate evidence for modification of transfer function parameters 
or definitions. 

8. Process the new set of measurements to estimate the original signals 
(prior to mixing) and store results of these processes as necessary. 

9. Evaluate successful separation and discrimination indicators. 

1 0. If memory allocation is insufficient or inappropriate go to step 1 . 

11. If modification of the filter formulation is adequate continue, else go to 
step 2. 

12. Go to step 3. 

One can translate the transfer function onto a state space representation with the added 
benefit of including the effect of initial conditions. This inclusion becomes especially 
important when one considers the mixing environment model in the general sense, 
e.g. the actual characteristics of the mixing transfer functions can vary in time. This 
particular method of the invention is described in the next section. 

State space time domain architecture and method 

FIGURES 6A and 6B show two mixing models for the state space time 
domain architecture. FIGURE 6A shows a general framework. FIGURE 6B shows a 
special case where A and B and are fixed, and its relation to conventional signal 
processing. Both models apply to multiple types of separation architectures. 

FIGURES 7A and 7B show two signal separation models for the state space 
time domain architecture. FIGURE 7A shows a general model and architecture. 
FIGURE 7B shows a special case, only the model is shown without the arrows in 
FIGURE 6A which depict parameter update procedures. 

This architecture of the invention models both the mixing and the separation 
environments as a linear dynamic system. In this way the parameters of the mixing 
environment are represented by the realization p = (A f B, C, D) as depicted in 
FIGURE 6A; whereas, p = (A, B, C 9 D) is the realization of the separation 
environment as depicted in FIGURE 7A. These parameters dictate the dynamics of 
the environments for mixing and separation. While the realization p is quasi- 
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constant and models the behavior of the mixing environment the realization p is to 
be defined by the separation process. 

The mixing model for the state space time domain architecture applies to 
multiple types of separation architectures. The model is an architectural depiction of 
the mathematical equations: 



Ax* Mm Equation (18) 



m = Cx + Ds Equation (19) 

On the other hand, as shown in FIGURE 7 A, the state space time domain 
signal separation architecture can mathematically be described as: 

x = Ax + Em Equation (20) 



u = Cx + Dm Equation (21) 

This architecture defines the necessary procedure for the utilization of this 
mathematical structure in signal separation, discrimination and recovery. This 
requires the definition of a procedure or a set of procedures for the signal separation, 
discrimination and recovery algorithms in such a way that when these procedures are 
exercised, the necessary parameters converge to some stable or quasi-stable solution, 
the outputs of the separation and discrimination process are the replicas of, or 
sufficiently similar copies of, the original signal sources S(t) . 



The method of this architecture of the invention: 

1 . is capable of processing of multiple, i.e. more than two mixtures 

2. is capable of including the effect of initial conditions directly 

3. may be extended to include time variance by assuming that the 
matrices (A, B 9 C, D) and (A, B, C, £>) are functions of time 

4. may allow for incorporation of general nonlinear models, e.g. for the 
mixing 
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x-Ifc-.IV Equation (22) 



m = $(x, 5, P 2 ) Equation (23) 

and for the separation 
x = F(x, W x ) Equation (24) 



u = <J>(x, m, ) 



Equation (25) 



5. may be extended to include nonlinear and time variant models by 
assuming that the operators I\ f , $ in Equations 22-25 are 
functions of time. 



In the linear case, one can show the inversion relationship between 
(A, B, C, 5) and (A, B 9 C, D) as follows: 
By Equation (19) 

s = -D ~*Cx + D ~ l m Equation (26) 

is the estimated source 

Equation (27) 



and therefore with reference to Equation (21) given that u 
signals s , 

A 

C = -D' l C 



and 

A 

D = -D 1 

Similarly, by substituting * = D~ x Cx + 
that 



Equation (28) 



1 w into Equation ( 1 8), one can show 
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A-A-BD-'C Equation^) 



and 

B tsD~ l Equation (30) 



In the event that D is not a square matrix, and is thus not invertible, such as the case 
where more measurements or mixtures are obtained than sources, one can follow a 
generalized procedure using the pseudoin verse [D T D]~ l assuming that [^D] , 
which is a square matrix, is in fact invertible. 

In this fashion, one needs to rebuild the network with the appropriate update 
laws so that the network parameters become the following alternate versions of 
Equations (26-30): 

s = -[5 T D] l D T Cx + [E^5y l D T m. Equation(31) 



^ -r - r - Equation (32) 

C = -[u D]D T C n 



D t-^DJD' Equation(33) 



- -t - i -t - Equation (34) 



r - i -t Equation (35) 



One may solve for the parameters of A , B , C , and Z> to attain their solutions 
precisely or approximately by a variety of suitable computing techniques and 
appliances, including a variety of software and hardware components endowed with 
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the practice of the invention. The input to such a computing appliance is an 
appropriate set of mixtures received through a plurality of detectors, with each 
detector receiving signals from a plurality of sources. After the processing of the 
input with the practice of the invention, the output from the computing device, is a set 
of separated signals, including restored source signals which closely estimate the 
original signals without interference. 

Using several families of information maximization criteria, a variety of 
optimization functions and numerical solution methods, both continuous (as in analog 
electronics) and discrete time (as in digital computing) may be used as appropriate, 
such that the signal separation problem is solved by computing the selected set of 
optimization functions for a given set of mixtures. 

One procedure for the utilization of this architecture for signal separation 
discrimination and recovery may be summarized as follows: 

1 . Construct and allocate a memory structure which is used for storing the 
measurements as well as the results of the data obtained from 
processing of the said measurements, and set measurement data 
acquisition or sampling rates. 

2. Construct four matrices, A , B , C , and D , of appropriate 
dimensions. 

3. Receive a new set of measurements, or mixed signals. 

4. Process the incoming mixture signals and measurements by performing 
signal conditioning, storage and related procedures. 

5. Evaluate signal separation criteria and functions. 

6. Update the elements of the four matrices, or a subset of them, 
accordingly. 

7. Accumulate evidence for modification of matrix dimensions, 
convergence rates, time constants and other parameters or definitions. 

8. Process the new set of measurements to estimate the original signals 
(prior to mixing) and store results of these processes as necessary. 

9. Evaluate successful separation and discrimination indicators. 

10. If memory allocation or data measurement sampling rate is insufficient 
or inappropriate go to step 1 . 
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11. If modification of the filter formulation is adequate continue, else go to 
step 2. 

12. Go to step 3. 

Special case. One can design a special case, where the canonical form of the 
mixing and separation model are modified from their canonical forms in Equations 
(18-21) so that A, A, B 9 B are fixed. In FIGURES 6B and 7B, this special case is 
illustrated pictorially. 

In conventional signal processing terms and related techniques, in relation to 
the mixing model, this special case yields a set of filtered and/or delayed versions of 
the source signals (s), upon which the mixing matrix C operates. The result of this 
operation is added to the result of the operation Ds . Thus, one can rewrite Equation 
(10) as 

M(s) = H(s)S(s) = C (si - AY l BS(s)B + DS(s) (Equation 36) 

Similarly, in relation to the separation model and architecture, this special case yields 
a set of filtered and/or delayed versions of the mixtures m , upon which the matrix 

C operates. The result of this operation is summed with the result of the 
operation D on m to obtain the separated signals u . 

U(s) = H(s)M(s) = [C{sI-AY x B+D]M(s) (Equation 37) 

The special case explicitly brings the congruence of the state space time domain 
architectures with the transfer function models. The procedure for the utilization of 
the canonical architecture form is very similar; the differences are highlighted as 
follows: 

1. In step 2, A and B can be fixed. 

2. In step 6, only the two matrices C and D would need to be updated. 
Ma pping architecture and method 

The mapping architecture and method is an alternative embodiment for 
facilitating the physical realization and implementation of the previously described 
architectures and methods, especially by hardware accelerated emulation or a digital 
signal processor device. Towards this goal, the above models can be viewed in 
discrete time form by simply replacing the time derivatives by samples in time, 
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i.e. r(/) - z(f+l) » where z can be x ox x in the above formulations. Below a 
complete procedure is outlined for the implementation of the mapping architecture 
and method: 

Let the mixing environment be represented by the time domain relation: 

L _ N 

m(t) = Ds(t) + £<V ( *>(0 - Y, A k mik) W (Equation 38) 

where (k) represents the kth derivative in the continuous time case, k delay samples in 
the discrete time case. As before, s(t) is the source, m(t) is the mixture, and u(t) 
are the separated output signals. 

The signal separation network will be defined by 

l' N* 

u(t) = Dm(t) + £c t m (4) (r) - J^A k u^(f) (Equation 39) 

where the dimensions of C and D are n x m, and the dimensions of A is n x n. 
Generally, L'^ L and N ? ^ N, but the case where L 1 < L and N f < N may also be 
considered. Thus the realization can now be expressed as: 

u(t) = W® (Equation 40) 

where 

W = [DC,...C L ,; -A r ..-A M i] (Equation 41) 

and 

«> r = [m(tfm < l > r ...m ^ u< 1)T ...ti^'] (Equation 42) 

Note that W is not necessarily a square matrix, i.e. W is of dimension n x [L'm + 
N f n]. 

The mapping model corresponds to the IIR (infinite impulse response) models 
in signal processing. Moreover, the special case of FIR (finite impulse response) 
models are obtained by eliminating the terms associated with derivatives of the output 
signal vector u{t) . 
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Physical realization and implementation 

Software emulation, the architectural representations of the transfer function- 
based frequency domain architecture, the state space time domain architecture and the 
mapping architecture may readily be translated into a computer program. For 
example, procedures above may be coded in a computer interpretable language and 
the resulting set of computer instructions may be executed on a compatible computing 
platform. Thus, software emulation is among the implementation options for this 
invention and involves the coding of the algorithm in a particular language in a way 
that it makes it possible to automatically execute it on a particular processor. High 
level languages such as FORTRAN and C/C++ have supporting processor-specific 
and general compilers that allow the code to become portable. Some of this code may 
be bundled with other functions and embedded onto an existing platform or included 
in processor's main program memory. On the other hand, more recent languages such 
as Java are platform independent so that the algorithm may be coded to run on 
multiple platforms in a network environment. 

Hardware accelerated emulation. The procedures may also be coded or 
translated from a higher level language (e.g. C/C++, FORTRAN, etc.) to a specific 
signal processing device or by embedding the procedure within a specific device 
program or memory device both inside or outside a signal processor. This constitutes 
hardware accelerated implementation of the architecture described. Hardware 
accelerated emulation implementation could be fast enough for real time on line signal 
separation, discrimination and recovery. The devices that may be used for this 
purpose include, but are not limited to emerging high power digital signal processors 
(DSP), microprocessors, application specific integrated circuits (ASIC), 
programmable devices, including but not limited to field programmable gate arrays 
(FPGA), reconfigurable devices, in system programmable devices, programmable 
logic arrays (PLA), or other custom built or programmed automated structures with 
pre-specified or custom processor cores that compute fast multiplication and addition 
and offer flexible high precision. 

FIGURE 8 shows one possible DSP implementation architecture 800. DSP 
architecture 800 includes one or more A/D (analog to digital) converters 802 
connected to a data stack 804. Data stack 804 is connected to a DSP 806 device, 



20 



WO 98/58450 



PCT/US98/13051 



which in turn has access to a memory device 808 and D/A (digital to analog) 
converters 810. The internals of DSP device 806 may include a variety of functional 
units as shown below. Different configurations are possible depending on the nature 
of the application, number of mixtures, desired accuracy, etc. , 
A block diagram depiction of a typical DSP device 806 is shown in FIGURE 
9. It may be convenient to integrate the signal separation process as part of other 
procedures already included to perform other signal processing functions particular to 
the device. A/D converters 802 and D/A converters 810 may be integrated into this 
structure for single chip solutions. DSP device 806 may include a memory device 
902, a MAC (multiplier accumulator ) 904, an ALU (arithmetic logic unit) 906, 
accumulators 908, shifters 910, internal accumulators 912, and address units and 
registers 914. Memory device 902 may include a program ROM data/program RAM. 
MAC 904 may include one or more multipliers and adders. ALU (arithmetic logic 
unit) 906 may include one or more adders, logic operators, and comparators. 
Accumulators 908 may include TEST circuits, serial ports, timers, wait state 
generators, clock generators, a host port, a host interface, stacks, and auxiliary 
registers. 

One typical procedure is as follows: The incoming signals are transduced 
electrically and these mixtures are digitized and stored. The procedures, algorithms 
and architectures of this invention are then carried out for the separation, 
discrimination, and recovery of the individual source signals. One or more of the 
units described, e.g. DSPs and/or DSPs with multiple functional units can be used to 
carry out this process for real time operations. The DSP units may be programmed to 
carry out the steps needed to solve the equations that are associated with the 
algorithm. This program may be obtained by compiling the high level language to 
translate the procedure into assembly or machine language which the DSP device can 
decode and carry out. Custom crafted machine or assembly language code may 
further accelerate the execution of the algorithm by optimizing the procedure for the 
DSP functional units. 

The outputs from these processes which form the approximations to the source 
signals may then be stored, converted to analog values, and/or processed further. 
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Examples of further processing include signal conditioning, classification, 
recognition, demodulation, and other operations for signal processing. 

Additional procedures that handle overdetermined (where number of mixtures 
are greater than the number of sources), underdetermined (where number of mixtures 
5 are less than the number of sources), varying (where the number of sources and 

mixtures vary) or unknown (where number of mixtures and/or sources are unknown) 
have been outlined in flowcharts. Furthermore, additional error checking, data 
integrity, and data identification schemes may enrich the possibilities and the 
application domains as well as robustness and reliability of the approach. 
10 These procedures may readily be programmed into the DSP or the device in 

which the necessary signal separation, discrimination, and recovery functions are 
embedded. 

Functions and crite ria for parameter adaptation 
Several families of adaptation techniques may be applied to obtain the 
1 5 parameters of the architectures and methods of this invention. 

Neurally inspired static architecture. As described by Equation (3) and 
FIGURE 2, signal separation in the case of a neurally inspired static architecture 
requires the computation of a weight matrix, W . For this, we first introduce the 
outer product matrix Q as 
20 ^ ^ v , - Equation (43) 

where / and g are two appropriate odd functions, which may or may not be 
related, and u is the set of output signals that estimate the input source signals u . 
This matrix is a function of time as u is a function of time. Further dependence of 
time can be built into the functions / and g . 
25 One can consider several categories of adaptation procedures based on Q in 

Equation (43), that yield the W matrix, by integrating one of the differential 
equations that mathematically describe the derivative of W , i.e. W . As an 
example, in discrete time digital implementation, W can be computed by Euler 
approximation that relates the values of W to W as: 

30 W t ^ = W t + hW t Equation (44) 
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where W t is the derivative of W with respect to time at time t, and h is the step in 

time between the two successive values of W , i.e. W x and W . 

Nine categories of weight update rules are described for W , as follows: 

W = T](a/-g) Equation (45) 

W = r\(aI-Q)W Equation (46) 

W = T\{al-Q)W T Equation (47) 

W = T)(adiag(Q)-Q) Equation (48) 

W = r\(adiag{Q)-Q)W Equation (49) 

W = r\(adiag(Q)-Q)W T Equation (50) 



W = madiag(Q)-Q) Equation (51) 



W = T)(arf/ag(e)-e)^ Equation (52) 



= Ti(a^(fi)-fi)^ Equation (53) 

where a is a positive number. In Equations (45-53), r\ is the rate of adaptation, where 
ti is a number that may vary during the process. The superscript (-T) represents 
inverse transpose, and diag(Q) is a diagonal matrix where all elements except the 
diagonal elements are equal to zero and the diagonal elements of diag(Q) are equal 



to those of Q . diag(Q) is time averaged values of the diag{Q) . 

State space architecture. As described by Equations (20-21), signal separation 
in the case of the state space time domain architecture requires the computation of 
matrices A , B , C , and D . For this, we first introduce two outer product 
matrices K and L as 

K =Au)g(x) T Equation (54) 
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L = fbQg(u) T Equation (55) 

where / and g are two appropriate odd functions, which may or may not be 
related, and u is the set of output signals that estimate the input source signals s , 

m is the set of mixtures received, and x is the set of internal states. These 
matrices are functions of time as u , m , and x are functions of time. Further 
dependence of time can be built into the functions / and g . 

One can consider several categories of adaptation procedures based on K 
and L in Equations (5 1-52), that yield the C and D matrices, assuming that, as 
given in the special case of the state space time domain representation, 
matrices A and B are fixed. As before with W , by integrating one of the 
differential equations that mathematically describe the derivative of C and D , one 
can obtain C and D . An example was given in Equation (44). 

Nine categories of weight update rules are described for C as follows: 

C = T\{yl-K) Equation (56) 

C = Ti(7 1-K)C Equation (57) 

C = I\(yl-K)C' T Equation (58) 

C = T\(ydiag{K)-K) Equation (59) 

C = tj(y diag(K) -K)C Equation (60) 

C = i\(ydiag(K)-K)C' T Equation (61) 



C = i\(ydiag(K)-K) Equation (62) 



C = v\(ydiag(K)-K)C Equation (63) 



C = *n(y diag(K) -K)C ~ T Equation (64) 
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where y £ 0 and diag{K) is time averaged values of the diag(K) , which may be 
realized by averaging one or more samples of the diag(K) matrix. 

In Equations (56-64), ti is the rate of adaptation, where r\ is a number that may 
vary during the process. The superscript (-T) represents inverse transpose, 
5 and diag(K) is a diagonal matrix where all elements except the diagonal elements 
are equal to zero and the diagonal elements of diag(K) are equal to those 
of K . K was described in Equation (54). 

Correspondingly, nine categories of weight update rules are described 
for D as follows: 

10 D = T)(a/-Z) Equation (65) 

D = r\(aI-L)D Equation (66) 

D = r\(aI-L)D' T Equation (67) 

D = r\(adiag(L)-L) Equation (68) 

D = r\(adiag(L)-L)D Equation (69) 

15 D = r\(adiag(L)-L)D T Equation (70) 



D = r\(adiag{L)-L) Equation (71) 



D = y\{adiag{L)-L)D Equation (72) 



D = r\(adiag(L)-L)D' T Equation (73) 



where a > 0 and diag(L) is time averaged values of the diag(L) , which may be 
20 realized by averaging one or more samples of the diag[L) matrix. 

In Equations (65-73), T| is the rate of adaptation, where ti is a number that may 
vary during the process. The superscript (-T) represents inverse transpose, 
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and diag(L) is a diagonal matrix where all elements except the diagonal elements are 
equal to zero and the diagonal elements of diag(L) are equal to those of L . L was 
described in Equation (55). 

Equations (56-64) and (65-73) could be combined in a variety of ways. Since 
nine categories are shown for C and D ,81 permutations are possible. 

Mapping architecture. The equations that apply to the neurally inspired static 
architecture, i.e. Equations (45-53) also apply to the mapping architecture, the 
parameters of which were described as a matrix W in Equations (40-41). To keep 
the dimension notation consistent, we rewrite the nine categories of weight update 
rules of this invention for the mapping architecture case. First, redefine L as: 

L =ftu)g(&) T Equation (74) 

where the dimensions of flu) is n x 1, and the dimensions of g($) r is 1 x N, where 
N = (L'+l)m + N'n in the notation of Equation (42). 

Then, the following categories of weight update rules can be designed: 

W = ti([/]0] -L) Equation (75) 

W = t]([/|0] -L)W Equation (76) 



W = TK[/|0]-W r Equation (77) 

W = i)(diag([L\0]) -L) Equation (78) 

W = r)(diag([L\0])-LW Equation (79) 

W = Ti(rf/ag([I|0]) -L)W ~ T Equation (80) 



W = T)(diagaL\Q])-L) Equation (81) 



W = r)(diagai\0])-LW Equation (82) 
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W = r](diag([L\0])-L)W T Equation (83) 

In Equations (75-83), r\ is the rate of adaptation, where r\ is a number that may vary 
during the process. The superscript (-T) represents inverse transpose, and diag(L) is 
a diagonal matrix where all elements except the diagonal elements are equal to zero 
5 and the diagonal elements of diag(L) are equal to those of L . The notation [I/O] 
refers to a non-square matrix of dimension nxN composed of positive number alpha 
times the n x n identity matrix / and the rest of the columns (i.e. columns N-n to N) 
are filled with zeroes. The notation diag[IfO] signifies the diagonal elements of the 
n x n submatrix obtained by taking the first n columns from L . Bars on a quantity 

1 0 signify time-average over one or more samples. W ~ T signifies transpose of the 

pseudoinverse of the non-square matrix W . Observe that the update equations of 
this generalized W can be straightforwardly factored out in terms of the update laws 
for the components matrices D , C , C L ,and A x ... A N i . 
Real time on line parameter adaptation 

1 5 The architectures and procedures described above can be augmented with 

additional techniques to improve the signal separation and recovery procedures. 

The unmixing (or, adaptive network weight) matrix entries may be updated for 
single instants of arbitrary data points calculated or obtained. Although by this 
method of application, it is possible to achieve real time on line signal separation, the 

20 method is also highly prone to errors since it allows a single point to alter the 

trajectory of the evolving solution. Another procedure for updating the parameters of 
the separation process applies the criteria to the entire data set, or selected data points 
from the entire data set. This method lacks causality essential to practical 
implementation: The related adaptation process does not progress in time or per 

25 sample, but utilizes the whole data set, including values in superseding instants. 

Moreover, a constant, static mixing matrix is assumed to apply over the whole range. 
Although this method is somewhat more robust than the first, it is essentially an 
offline method not suitable for real time signal separation. Furthermore, when the 
assumption of a static constant matrix is incorrect, the accuracy of the unmixing 

30 process suffers. 
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Continuous multithreaded adaptation. FIGURE 10 shows multithreaded 
adaptation processes. Several signal separation processes may be overlapped in time 
to prevent parameter saturation and reduce the likelihood of sustained incorrect 
parameter estimation. Each rectangle represents a separation process. At the 
5 beginning of each process, the parameters are initialized and adaptation criteria are 

applied. At each point in time, there is one master process, the separation results from 
which are reported, and one or more slave processes. 

Thus at each point in time, there is more than one process, the separation 
results from which are reported. Depending upon the application, each process output 

10 may be used in a variety of ways. Suppose, for example, among separated signals, 
only a single signal of interest needs to be estimated. In that case, there would be a 
master process that optimizes the quality of that signal of interest, whereas other 
process results are ignored. In another example, the redundant (or slave) process 
results may be compared with those from the master process and/or those from other 

15 slave processes to (1) access the success of each separation process, (2) evaluate and 
select the best initialization points, or (3) decide when a process should expire. When 
a master process expires, a slave process takes its place. In FIGURE 10, at time index 
t„ process 1 expires. At that point, depending upon the quality measure of their 
separation results, the best one of the other three slave processes may take its place as 

20 the master process. 

Formation of a data buffer or set (long stack) to be used for parameter 
adaptation. As discussed, the application of signal separation criteria and computation 
and adaptation of separation parameters are executed on either a single data point or 
on the whole data set obtained. Both of these are of limited use. In the case of the 

25 single data point, it is possible and probable that a single or a series of noisy points 
would change or divert the trajectory of the adaptation parameters from their ideal 
solution. On the opposite end of the spectrum, applying the separation criteria to the 
whole data set is often impractical since (i) the storage of every data point may not be 
possible and (ii) the latency created by the processing of all data points may be 

30 unacceptable. 

For this reason, the present invention involves the creation of a real time on 
line running time window the size and contents of which can be adjusted depending 
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upon various criteria imposed upon the measured signals, separated signals, or the 
time constraints of the particular application. These criteria would generally measure 
the level of success of the separation or discrimination process, such parameters as 
power spectrum of the measured or separated signals, statistical similarity of the 
separated or measured signals, or be guided by the human user of the separation 
process. The long stack may be contained in a storage device, a digital or analog 
recording medium. 

The contents of the data set may be modified after a new set of measured 
signals are obtained. The most straightforward way this may be done is by simply 
eliminating from the data set the oldest data points. This may be done by constructing 
a first in first out (FIFO) type stack as shown in FIGURE 1 1 . FIGURE 1 1 shows one 
embodiment of a constant width long stack based on a first in first out (FIFO) stack 
construct and data acquisition process. The shaded region indicates the contents of 
the long stack in this case. Until sample time = W, every data set obtained is stored. 
When sample time = W, the stack is full. Beyond this point, the oldest data set is 
deleted from the stack. Note that this is a simplified illustration of the invention and 
the stack may not be of FIFO type or of constant width. The procedure below is an 
example of a program based on this invention for selecting a data set: 

1 . Set the length of data stack. 

2. Initialize parameters of signal separation. 

3. Obtain the set of measurements of the current time instance. 

4. Is this a usefiil data set or point? 

No -> Is it desirable to perform separation on this data set? 

Yes — ► Go to step 6. 

No-* Go to step 3. 
Yes -+ Is the data stack full ? 

No — ► Store data point as the most recent data point. 

Yes — ► Discard oldest data point store this data point as the 
latest. 

5. Is it desirable to perform separation on this data set? 
Yes -* Go to step 6. 

No-> Go to step 3. 
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6. Perform the signal separation procedure. 

7. If this was a useful data point then store the results of the separation 
process in the data stack. 

8. If checkpoint is reached, evaluate the success of the separation process. 
5 9. Is separation satisfactory? 

No-* Go to step 1. 

Yes — ► Is it desirable to reduce stack size? 
No-+ Go to step 3. 
Yes-* Go to step 1. 

10 It is possible to introduce latencies into the process, that is one may output the 

result of the signal separation performed on samples obtained prior to the current time 
instance, rather than the sample most recently available. 

Stochastic or deterministic selection of a subset of elements (short stack) from 
the long stack for application of adaptive parameter estimation criteria. Prior to the 

15 execution of the signal separation and discrimination algorithms, in particular, prior to 

the application of the adaptive parameter estimation criteria, one needs to select which 
data points from the long stack to use. A simplified illustration of the concept is 
shown in FIGURE 12. FIGURE 12 shows an example of a constant width long stack 
based on a first in first out (FIFO) stack construct that is seven elements deep. The 

20 shaded region indicates the contents of the long stack in this case. Until sample time 
= 7, every data set obtained is stored. When sample time index = 7, the stack is full. 
Beyond this point, the oldest data set is deleted from the stack. Two elements are 
selected from the long stack to form the short stack at each process time index. These 
elements are shown as black squares. 

25 No such selection is required if one uses either the latest set of results and 

measurements or the entire set of results and measurements. The latter computation is 
noncausal since present outputs will depend on future mixtures. Hence the latter 
computation can only be performed offline. As already presented in real time 
execution of the signal separation and discrimination algorithms this is often not 

30 practical. 

It may be desired that one uses a select set of the entire data stack available - 
either by deterministically, randomly, or pseudorandomly selecting a set of points 
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from it. The deterministic selection is entirely predictable and will be the same with a 
set of identical data each time. The use of random or pseudorandom selection criteria; 
however, may yield different selected sets even with identical data. 

One powerful deterministic criteria is to expand or contract the size of the 
5 short stack starting from, expanding around or ending at a random or pre-determined 
index of the long stack. This process may be guided by the results of the separation 
process, i.e. the quality of the separation and discrimination received may be used for 
defining the index of the data in the long stack, increasing or decreasing the number of 
data elements of each point in the short stack, the size of the short stack and the 

10 numerical accuracy of the data in the short stack. Deterministic procedures may be 
combined with stochastic procedures, such as those that randomly select the length, 
numerical accuracy, starting point and elements of the short stack. The stochastic 
procedures may utilize appropriate probability distribution functions to define the 
indices of the elements in the long stack that will be included in the short stack. For 

15 instance, the probability distribution of indices may favor or disfavor the selection of 

recent points. The definition of the probability distribution may itself be a stochastic 
or deterministic process based on original and separated signal and mixing 
characteristics at the time of the computation. Computations of updates for the 
network may all be causal, i.e. the present value of the output depends only on present 

20 and past values of the mixture. 

FIGURE 13 shows the application of the separation criteria using a set of data 
points (short stack) from a long stack which contains all relevant sets of data points, 
measured or computed over a window of time. In FIGURE 13, it is shown pictorially 
the construction and elements of the long and short stacks. The long stack is N 

25 elements wide and M elements deep. The functions are grouped by their time index k. 
There are N time instances which contain M results of the computation, including the 
measured signals and the estimated original sources or the separated signals, the 
results of intermediate computations, and various functions thereof. The contents of 
the long stack represent a moving window or history of measurements and 

30 computations. The most recently measured and computed data points are appended to 
the stack entry with the latest index whereas the least recent set of said points are 
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purged from the stack. As suggested earlier, the size of the stack itself need not be 
constant. 

The short stack is a subset of selected entries from the long stack for the 
computation at a particular time instance. The short stack contents are also modified 
5 at every instance in time when a computation is performed. The short stack entries 
selected at time t are used at time t to compute the estimated original signals or 
separated signals. However, note that although time t is the time of computation, it 
need not be the time associated with the separated output. There may be a latency in 
the output that is a look-ahead scheme could be at work so that the computation at 
10 time t is the separated signal for time t - L where L represents the latency. 

Thus, the algorithm of the present invention and framework consider general 
time domain systems making the invention capable of accommodating changing 
phenomena of signal delays and other nonlinear events which are routinely 
experienced in practice with signal and wave propagation, transmission and 
15 transduction. 

The present invention also allows for statistically significant intervals to be 
processed and exploits the various forms of weighing of the outputs by manipulating 
the criteria of selecting the entries of the long and particularly the short stack. 
Recursive procedures for adaptive parameter estimation. In a variety of 
20 applications, especially when high accuracy of the separation process is desired, one 
. may recursively apply the signal separation process, that is the outputs from the signal 
separation process can be treated as inputs to the next separation process until desired 
accuracy is achieved. Each successive process of separation may use the same or 
different model, architecture, and parameters from the previous process. This method 
25 may also be combined with other real time on line procedures of this invention, such 
as the multithreaded adaptation procedure. 
Audio signal processing application 

For speech, speaker, and/or language recognition and processing, given the 
computational challenges associated with these processes, or due to objectives of user 
30 privacy and/or surveillance, it is desirable to obtain noise free clear signals, e.g. 

speech segments that are uncorrupted by noise, environmental sounds, and other 
speakers. The signal separation architectures and algorithms of this invention may be 
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used as a front end towards this goal. Microphone arrays may be used to obtain 
multiple versions of speech (and other) signals. These may be the mixture signals 
used in signal separation to obtain the original individual signals or speech segments. 
These segments will contain less noise and other sources of interference than the 
original signals picked up by the microphones. 

FIGURE 14 shows an audio application based on the signal separation and 
recovery procedures of this invention. Audio signals are converted electrical signals 
by the elements of a microphone array 1402. Each element of microphone array 1402 
receives a different version (or mixture) of the sounds in the environment. Different 
arrangements of microphone elements may be designed depending on the nature of the 
application, number of mixtures, desired accuracy, and other relevant criteria. 
Following some signal conditioning and filtering, these mixture signals are converted 
from analog format to digital format, so that they can be stored and processed. The 
DSP device 806 of the system is programmed in accordance with the procedures for 
signal separation and recovery procedures of this invention. The internals of DSP 
device 806 may include a variety of functional units for various arithmetic and logic 
operations, and digital representation, data storage and retrieval means to achieve 
optimum performance. Circuits and structures shown in figure may undergo further 
integration towards realization of the whole system on a single chip. FIGURE 14 
illustrates the application specific system in the framework of the hardware 
accelerated implementation option shown in FIGURE 8. 

FIGURE 15 shows three audio input interfaces to the DSP based handheld 
device that may be embedded into the microphone-DSP interface. On the very left is 
the prevalent interface to a microphone 1502. Microphone 1502 is connected to an 
A/D converter 1504, which is in turn connected to a DMA (direct memory access) 
channel 1506 and DSP or voice processor chip 1508. In the middle it is shown that 
microphones 1 502 in a straightforward way, as is the case in FIGURE 14. This may 
use another A/D converter 1504 and perhaps more importantly another DMA channel 
1506 of the DSP or processor chip 1508. On the very right we show a potential way 
in which this impact could be reduced using a multiplexer 1510. The select signal to 
multiplexer 1 510 is not shown, but it may be generated from the sampling clock in a 
straightforward way. Note that in the case of the multiplexed input, the impact on the 
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audio signal path is minimized. One may also consider a slight variation where, if 
simultaneous sampling is desired, one may sample using two A/D converters 1504 
and multiplex the digitized data instead of the (analog) microphone outputs. 

FIGURE 16 shows a schematic of audio device suitable for use with the 

5 present invention. This is a smart sound sensing and processing device 1600 that may 
be programmed to perform a variety of functions. This device may be used in whole 
or in part within a computing or communications device to implement voice or 
acoustic signal based interface. Smart microphone 1600 receives a plurality of sounds 
from sources 1602 through a microphone array 1604. Microphone array 1604 is made 

10 up of an array of microphone elements 1605. Microphone array 1604 is connected to 
an A/D converter 1606 and one or more DSPs or integrated DSP cores 1608. DSP 
core 1608 is connected to a memory device 1610, at least one digital circuit 1612 and 
analog circuit 1614, an online programming interface 1616, and a D/A converter 
1618. Memory device 1610 may contain the signal separation and recovery 

15 algorithms and procedures described earlier. Digital circuit 1612 and analog circuit 
1614 may include signal conversion, conditioning and interface circuits. D/A 
converter 1618 feeds output channels 1 620 and other outputs 1 622. Although it is 
shown that they lay on a flat surface in the drawing, the microphone axes need not be 
coplanar. In fact, a cylindrical arrangement similar to the one shown in FIGURE 17 

20 would be desirable to pick up sounds from all directions. Some of these directions 
may be disabled. This smart microphone contains microphone elements 1605 
depicted in its schematic in FIGURE 16. The device is minimally capable of audio 
signal separation and recovery tasks. Additional features may be added by 
programming these features as well onto the digital signal processor. This device may 

25 be used as a replacement for ordinary microphones after being pre-programmed or 
used in new audio interfaces that exploit its flexibility by online programming. 
Although it is drawn larger than the conventional microphone to illustrate its 
components, it may actually be of comparable size or smaller due to its high level of 
integration. 

30 Microphone technology and microphone arrays. A microphone is a device that 

converts acoustic energy received as vibratory motion of air particles into electrical 
energy sent along the microphone cable as a vibratory motion of electrons. Once this 
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conversion has taken place, the sound information is freed from the normal acoustic 
constraints. It can be amplified, transmitted on a wire or via radio waves, it can be 
stored, and processed. 

Today the microphone is a ubiquitous device. Many types of microphones 
exist and this section presents a brief overview of the technology associated with 
microphones. The desirable features of any microphone depend to some extent on the 
particular applications for which it is intended. However, in studying the principles of 
acoustic transduction as well as the manufacturer's specifications of microphones, the 
following performance categories have an important bearing on the selection of . 
appropriate microphones: 

1 . frequency response on axis, which refers to signal output versus 
frequency for a constant acoustic level 

2. directivity, which refers to the response of the microphone to sounds 
received from all angles on a plane, e.g. omnidirectional, cardioid, 
shotgun 

3. frequency response off-axis, same as (1) but off the off-axis 

4. sensitivity, or the efficiency of the microphone to convert acoustic 
energy to electrical signal 

5. self noise, or the inherent noise level 

6. distortion, in other words, the deviation from a linear response to 
acoustic energy 

. With the exception of a few exotic types based on heated wire or a cloud of 
ions, all practical microphones make the conversion of acoustic energy to mechanical 
energy via the mechanical vibrations of in response to sound waves of a thin, light 
diaphragm. Generally this diaphragm is circular in shape and clamped at its 
periphery, although other shapes do occur and including the thin ribbon variety 
stretched between clamps at each end. It appears therefore that the energy conversion 
in microphones commonly takes place in two stages, albeit simultaneously: acoustical 
to mechanical and mechanical to electrical. 

In the first stage, the two main ways in which a microphone of any transducer 
type extracts mechanical energy from the sound wave are pressure operation and 
pressure gradient operation. The distinguishing feature of pressure operated 
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microphones is that the rear surface of the microphone is enclosed so that the 
actuating force is that of the instantaneous air pressure at the front. A small vent is cut 
through the casing to equalize the long term internal and external air pressures. A 
purely pressure operated microphone is omnidirectional. A pressure-gradient 
microphone, on the other hand, is built with both faces of its diaphragm equally open 
to air. The effective force on the diaphragm at any instant is not simply due to the 
pressure at the front end but to the difference in pressure, or the pressure gradient 
between the front and the back. This has an important bearing on the system's 
directivity. A whole family of directivity patterns become possible when one 
combines pressure and pressure gradient operation. 

The second stage may use a variety of electrical generator principles to convert 
the extracted mechanical energy to electrical energy and microphones tend to be 
categorized accordingly. Some of the common categories can be delineated as 
follows: 

1 . moving coil (dynamic) microphone based on the principle that motion 
of a conductor in an electric field generates an EMF causing flow of 
current in the conductor 

2. ribbon microphone based on the same principle as the dynamic 
microphone, where the ribbon acts as the diaphragm and the conductor 

3. condenser (capacitor or electrostatic) microphone based where the 
transducer element is a capacitor which varies based on the vibration of 
the diaphragm 

4. electret microphone which uses a polarized or neutral diaphragm with 
the fixed plate coated with electret material (back polarized) 

5 . piezoelectric microphone using crystalline or ceramic materials which 
possess piezoelectric properties where vibrations of a diaphragm are 
transmitted to the bimorph (oppositely polarized slabs joined to form a 
single unit) by a rod giving rise to an alternating voltage at the output 
terminals proportional to the effective displacement 

6. the carbon microphone, used in many telephones, which uses granules 
of carbonized hard coal that when subjected to external pressure 
variations will cause their areas of contact to increase or decrease 
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During the past decade silicon micromachining techniques have been 
successfully applied to the fabrication of miniature microphones on silicon wafers. 
Silicon microphones based on different principles, e.g. piezoelectric, piezoresistive, 
and capacitive principles are also available. There is growing interest in this area due 
to numerous advantages in improved dimension control, extreme miniaturization, the 
ability to integrate on-chip circuitry, and potential low cost as a result of batch 
processing. 

Here we include discussion of two types of micromachined microphones, 
namely capacitive and piezoelectric. 

Micromachined capacitive microphones. Most silicon microphones are based 
on capacitive principles because of the higher sensitivity flat frequency response and 
low noise level. The capacitive microphones can be divided into electret 
microphones, condenser microphones and condenser microphones with integrated 
field effect transistors (FETs). The capacitive microphone consists of a thin flexible 
diaphragm and a rigid backplate. The two parts can be realized either as two separate 
silicon chips or one single chip. On the backplate acoustic holes are formed by use of 
anisotropic etching. The diaphragms of the silicon microphones can be thin Mylar 
foil, polyester, a low pressure chemical vapor deposition (LPCVP) silicon nitride film, 
etc. For the two chip structural microphone silicon direct bonding, anodic bonding, or 
polymer adhesives are used to assemble the diaphragm and the backplate. Such 
assembly progress always involves laborious alignment procedures, and given that 
most bonding processes also involve procedures that affect integrated electronics and 
change the material characteristics it seems that microphones with single chip 
solutions should be preferable. 

The fabrication process sequence of this microphone uses seven masks 
including silicon anisotropic etching and sacrificial layer etching. It is claimed that 
diaphragms with corrugations have several advantages for the design of silicon 
condenser microphones. First, the corrugation reduces the intrinsic stress of the thin 
diaphragm which results in higher sensitivity. Second, it can be fabricated by use of a 
single wafer process. Acoustic holes are produced automatically after the air gap is 
formed, i.e. the sacrificial layer is removed and the diaphragm is released. The 



37 



WO 98/58450 



PCT/US98/13051 



sensitivity and the resonance frequency of the diaphragm can be optimized by 
choosing appropriate structure parameters. 

Micromachined piezoelectric microphones. The performance limits of 
micromachined piezoelectric microphones has been explored both theoretically and 
experimentally. In addition, such devices have been manufactured and tested using 
residual-stress compensation with on-chip, large-scale-integrated (LSI) CMOS circuits 
in a joint, interactive processes between commercial CMOS foundries and university 
micromachining facilities. One such 2500 x 2500 x 3.5 jim cube microphone has a 
piezoelectric ZnO layer on a supporting low-pressure chemical-vapor-deposited 
(LPCVD), silicon-rich, silicon nitride layer. The packaged microphone has a resonant 
frequency of 18 kHz, a quality-factor Q = 40 (approx.), and an unamplified sensitivity 
of 0.92 fiV/Pa, which agrees well with the calculated sensitivity. Differential 
amplifiers provide 49 dB gain with 13 ^tV A-weighted noise at the input. 

Related processing techniques to align features on the front side of a wafer to 
those on its backside have been attempted for bulk micromachining. Towards this 
end, a tiny (30 jim-square and 1.6 |nm-thick) diaphragm can serve as an alignment 
pattern. At the same time that the alignment diaphragm is made, much thicker, large- 
area diaphragms can be partially etched using "mesh" masking patterns in these areas. 
The mesh-masking technique exploits the etch-rate differences between (100) and 
(111) planes to control the depths reached by etch pits in selected areas. The large, 
partially etched diaphragms (2 to 3 mm squares) are sufficiently robust to survive 
subsequent IC-processing steps in a silicon-foundry environment. This reported size 
will decrease with the availability of smaller feature sizes. The thin alignment 
diaphragm can be processed through these steps because of its very small area. The 
partially etched diaphragms can be reduced to useful thicknesses in a final etch step 
after the circuits have been fabricated. This technique has been successfully employed 
to fabricate microphones and on-chip CMOS circuits. 

Use of micromachined piezoelectric microphones creates the opportunity to 
build electrothermally tunable resonant frequencies. Specific devices with this 
property achieve resonant-frequency modulation, f m via the thermal expansion induced 
by on-diaphragm, polysilicon heaters. A near linear decrease in f m with increasing 
resistor power at the rate of -127 Hz/mW has been measured between 20.8 and 15.1 
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kHz. The mechanical quality-factors Q m are approximately 100. Moreover, the 
temperature-power curve is linear with slope 0.3 Celsius/mW at the hottest sense 
element on the diaphragm. The variable resonant frequency and high quality-factor 
provide an acoustic filtering capability which may have applications to ultrasonic 
range finders, velocity sensors, signaling elements, and speech processing. 

Unlike its predecessor with diaphragms clamped at four ends, combined 
microphone-microspeaker devices contain a cantilever that is free from the residual 
stresses found with clamped diaphragms. Use of the cantilever allows for sensitivity 
levels greater than other microphones with micromachined diaphragms. In addition, 
when the device is driven electrically as an output transducer, i.e. microspeaker, the 
relatively large deflections of the free end produce significant acoustic output. 

Theory of the integrated microphone can be developed through combining 
mechanics, piezoelectricity and circuit theory. Also required are theoretical 
optimizations for sensitivity-bandwidth product and signal-to-noise ratio. 

Microphone array assembly and adaptive tuning of microphone elements. The 
current use of the term "microphone array" refers to arrangements of microphones 
across large distances, in the order of meters. This is a result of two factors. First, the 
diameter of most common microphones do not allow for smaller array dimensions. 
Second, so far there has been little need to pack microphones very closely since the 
primary motivation of building microphone arrays is to provide the necessary input to 
beamforming algorithms, which, due to the large wavelength associated with most 
audio signals (0.1 - 1.0 meter) require large separation of microphone elements. For 
example, a one dimensional microphone array used for identifying speaker location 
includes 5 1 microphones are arranged in a line with a uniform separation of 4 cm 
across a distance over 2 m. This arrangement provides selectivity without spatial 
aliasing for frequencies up to 4000 Hz. 

Only with the advent of technology and new algorithms, with the 
understanding of the workings of the inner ear, we now have a reason to build 
structures resembling the inner ear, to do the kinds of signal processing needed for 
truly smart audio interfaces. From this point forward in this description, the term 
microphone array is used to refer to compact structures, less than 5 cm in any 
dimension. At this size, a single device or interface can contain the array and there is 
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no need to wire a whole room or enclosure with microphones in order to carry out 
sophisticated audio signal functions, such as signal separation. Although the physical 
dimensions make beamforming and localization mathematically impractical at this 
size, the ultimate goal is to integrate all the functions that a person can perform with 
one good ear onto such a device. Multiple devices can then be arranged to solve 
further problems that require physical dimensions suitable for source localization and 
beamforming. 

The integration of the microphone onto silicon enables tight control of 
microphone dimensions and spacing. This in turn is expected to result in more 
reliable exploitation of beamforming algorithms. However, the response 
characteristics of the micromachined microphone elements are less predictable than 
their conventional counterparts. Several issues involved in tight coupling of the 
piezoelectric microphones remain unresolved. Note that we are viewing the 
individual micromachined or microelectromechanical system (MEMS) device in this 
context as a single component in an array of similar units, analogous to a single pixel 
in a CCD array or a transistor on an integrated circuit. 

Since the signals from each MEMS device commonly exhibits offsets and 
nonlinearities due to existing physical constraints in the manufacturing process, the 
transduced signals must rely on off-the-shelf microelectronics for conditioning, 
testing, and one-time factory calibration in order to generate reliable measurements for 
diagnosis and control. Such one-time factory calibration, however, is not sufficient 
for generalized use due to, e.g. radiation, device aging, and temperature variations. 
Furthermore, when one considers tightly coupling similar and even disparate MEMS 
devices on the same substrate, external fixes to the problem would overcome the 
advantages of miniaturization afforded by MEMS. 

This application of the invention uses integrated adaptive systems alongside 
the acoustic sensors because it provides the ability to modify the relationships on or 
offline in the presence of anomalies and nonlinearity. In addition, using modular 
design structure renders the augmentation of the signal conditioning system 
straightforward. 

Digital signal processors. There exist many digital signal processors (DSP) for 
audio, speech, and multimedia applications. It will be up to the tradeoffs in a 
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particular application to select the best processor families for the separation and 
recovery of audio signals. The choices are between various architectures and 
numerical representations, e.g. floating or fixed point. To achieve a highly integrated 
solution (e.g. one chip) may require the embedding of a DSP core either from a pre- 
designed device or designed from standard silicon cell libraries. 

The compiler front-end to the DSP assembler and linker creates portability 
between the two implementation options described in this invention. There is thus a 
two-way connection between the high level language algorithm of the software 
emulation implementation and the hardware accelerated emulation implementation. 
In addition, a similar direct link exists between many computing environments and the 
DSP emulation environments, for example, C/C++ library and compilers for various 
processors. 

It is desirable to generate fast and compact assembly code specifically for each 
digital signal processor. Unfortunately, the assembly code generated by a specific 
DSP compiler is often not as optimal as desired. Various tool environments for digital 
processor design and co-design are composed of a hierarchy of models. Starting from 
the higher behavioral level, this layered description facilitates the identification and 
correction of design errors before they are physically committed. Moreover, it also 
gives a fairly accurate performance estimate provided that the described devices can in 
fact be physically realized. 

Programmable logic can be an integral part of the related development process. 
A programmable DSP core (a DSP processor that is designed for integration into a 
custom chip) may be integrated with custom logic to differentiate a system and reduce 
system cost, space, and power consumption. 

Signal conversion, conditioning and interface circuits. Signal conversion, 
conditioning and interface circuits are readily available as discrete components, e.g. 
A/D and D/A converters, memory ICs, timers, op amps, filters, etc. For the sake of 
rapid prototyping and programmability it may be desirable to use programmable 
devices, both analog and digital. 

Many digital (field programmable gate arrays) FPGA vendors are offering 
increasing complexity and gate and pin count devices for programming. Further 
reference is made to reconfigurable hardware in its various forms, which include 
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dynamic system resident reconfiguration and virtual hardware and reconfigurable DSP 
cores. For adaptive processing algorithms and their implementation, there exists a 
growing opportunity for evolving seamless integrations of software and hardware, as 
well as integration of processor and sensor. Given that the invention involves both 
analog and digital signals reference is also made to analog programmability of 
hardware, e.g. Electrically Programmable Analog Circuit (EPAC). It may be 
advantageous to combine analog, digital and mixed signal programmable devices with 
custom electronic models or discrete components as necessary. It may also be 
beneficial to pursue custom circuits and structures that can ultimately be integrated 
readily alongside acoustic sensing devices. 

The foregoing description of the invention has been presented for purposes of 
illustration and description. It is not intended to be exhaustive or to limit the 
invention to the precise forms disclosed. Obviously, many modifications and 
variations will be apparent to practitioners skilled in this art. It is intended that the 
scope of the invention be defined by the following claims and their equivalents. 

What is claimed is: 
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CLAIMS 

1 1 . A signal processing system for separating a plurality of input signals 

2 into a plurality of output signals, the input signals being composed of a mixture of a 

3 plurality of source signals that have been affected by a medium, the source signals 

4 being associated with a plurality of sources, the output signals estimating the source 

5 signals, the system comprising: 

6 a plurality of sensors for detecting the input signals; 

7 a storage device for receiving and storing the input signals; 

8 an architecture processor for defining and computing a signal separation 

9 architecture, the signal separation architecture defining a relationship between the 

10 input signals and the output signals, the relationship having constant parameters and 

1 1 time varying parameters; 

12 an update processor for computing a rate of change for each time varying 

13 parameter and the time varying parameters in response to the rate of change associated 

14 with each time varying parameter; and 

15 an output processor for computing the output signals based on the signal 

16 separation architecture, the constant parameters, and the time varying parameters. 

1 2. A signal processing system according to claim 1 , wherein the plurality 

2 of sensors are arranged in a sensor array, the sensor array having a directional 

3 response pattern. 

1 3. A signal processing system according to claim 2, wherein the 

2 directional response pattern of the sensor array is capable of being modified by 

3 performing signal processing on the input signals. 

1 4. A signal processing system according to claim 1 , wherein a quantity of 

2 the input signals and a quantity of the output signals are not equal. 

1 5. A signal processing system according to claim 4, wherein at least one 

2 output signal is a function of at least two source signals. 
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6. A signal processing system according to claim 4, wherein at least two 
output signals are functions of a same source signal. 

7. A signal processing system according to claim 1 , wherein the 
computing of the output signals is based also on a plurality of internal states of the 
system. 

8. A signal processing system according to claim 1, wherein the 
computing of the output signals is based also on at least one of the input signals, the 
output signals, previously received input signals, and previously computed output 
signals. 

9. A signal processing system according to any one of claims 1-8, 
wherein the signal separation architecture is defined by a state space representation 
that establishes the relationship between the input signals and the output signals. 

10. A signal processing system according to claim 9, wherein the 
computing of the output signals is based also on current states of the state space 
architecture. 

11. A signal processing system according to claim 9, wherein the 
computing of the output signals is based also on previously computed states of the 
state space architecture. 

12. A signal processing system according to claim 9, wherein the 
computing of the output signals is based also on current states and previously 
computed states of the state space architecture. 

13. A signal processing system according to any one of claims 9-12, 
wherein the state space representation is mapped onto a finite impulse response (FIR) 
filter. 
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1 14. A signal processing system according to any one of claims 9-12, 

2 wherein the state space representation is mapped onto an infinite impulse response 

3 (IIR) filter. 

1 1 5. A signal processing system according to any one of claims 9-12, 

2 wherein the state space representation is generalized to a nonlinear time variant 

3 function. 

1 1 6. A method for separating a plurality of input signals into a plurality of 

2 output signals, the input signals being composed of a mixture of a plurality of source 

3 signals that have been affected by a medium, the source signals being associated with 

4 a plurality of sources, the output signals estimating the source signals, the method 

5 comprising: 

6 receiving the input signals; 

7 storing the input signals; 

8 defining a signal separation architecture, the signal separation architecture 

9 defining a relationship between the input signals and the output signals, the 

10 relationship having constant parameters and time varying parameters; 

1 1 initializing the constant parameters and the time varying parameters; 

12 computing a rate of change for each time varying parameter; 

1 3 computing the time varying parameters in response to the rate of change 

14 associated with each time varying parameter; and 

15 computing the output signals based on the signal separation architecture, the 

16 constant parameters, and the time varying parameters. 

1 1 7. A method according to claim 16, wherein the rate of change of at least 

2 one set of time varying parameters is contained in an array having at least two 

3 dimensions. 

1 1 8. A method according to claim 1 7, wherein at least one of the arrays 

2 includes a function of the outer product of an expansive function of a set of the output 
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3 signals arranged in a first one dimensional array and a compressive function of a set of 

4 the output signals arranged in a second one dimensional array. 

1 1 9. A method according to claim 1 6, wherein the signal separation 

2 architecture is defined by a state space representation that establishes the relationship 

3 between the input signals and the output signals. 

1 20. A method according to claim 1 9, wherein the rate of change of at least 

2 one set of time varying parameters is contained in an array having at least two 

3 dimensions. 

1 21 . A method according to claim 20, wherein at least one of the arrays 

2 includes a function of an outer product of an expansive function of a set of the output 

3 signals arranged in a first one dimensional array and a compressive function of a set of 

4 internal states of the state space architecture arranged in a second one dimensional 

5 array. 

1 22. A method according to claim 1 7 or 20, wherein at least one of the 

2 arrays includes a function of an outer product of at least one expansive function of a 

3 set of the output signals arranged in a first array having at least two dimensions and at 

4 least one compressive function of a set of the output signals arranged in a second array 

5 having at least two dimensions. 

1 23. A method according to any one of claims 1 6-22, wherein a plurality of 

2 the methods are overlapped in time. 

1 24. A method according to claim 23, wherein at least one of the methods 

2 uses a predetermined set of constant values or a random set of numbers for initializing 

3 the parameters. 

1 25. A method according to claim 23, wherein at least one of the methods 

2 uses parameters computed previously by another method overlapping in time. 
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1 26. A method according to claim 23, wherein at least one of the methods is 

2 terminated. 

1 27. A method according to claim 26, wherein at least one method uses the 

2 parameters computed by previously terminated methods. 

1 28. An acoustic signal discrimination system for discriminating a plurality 

2 of signals into a plurality of output signals, the input signals being composed of a 

3 mixture of a plurality of source signals that have been affected by a medium, the 

4 source signals being associated with a plurality of sources, the output signals 

5 estimating the source signals, the system comprising: 

6 a plurality of acoustic sensors for detecting the input signals, the input signals 

7 being composed of a set of functions of a set of the source signals; 

8 a storage device for receiving and storing the input signals; 

9 an architecture processor for defining and computing an acoustic signal 

10 discrimination architecture, the architecture defining a relationship between the input 

1 1 signals and the output signals, the relationship having constant parameters and time 

12 varying parameters; 

13 an update processor for computing a rate of change for each time varying 

14 parameter and the time varying parameters in response to the rate of change associated 

1 5 with each time varying parameter; and 

16 an output processor for computing the output signals based on the acoustic 

17 signal discrimination architecture, the constant parameters, and the time varying 

18 parameters. 

1 29. An acoustic signal discrimination system according to claim 28, 

2 wherein the plurality of acoustic sensors are arranged in an acoustic sensor array, the 

3 acoustic sensor array having a directional response pattern. 

1 30. An acoustic signal discrimination system according to claim 29, 

2 wherein the directional response pattern of the acoustic sensor array is capable of 
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3 being modified by processing of the signals detected by the acoustic sensors of the 

4 acoustic sensor array. 

1 3 1 . An acoustic signal discrimination system according to claim 28, 

2 wherein a quantity of the input signals and a quantity of the output signals are not 

3 equal. 

1 32. An acoustic signal discrimination system according to claim 3 1 , 

2 wherein at least one output signal is a function of at least two source signals. 

1 33. An acoustic signal discrimination system according to claim 31, 

2 wherein at least two output signals are functions of the same source signal. 

1 34. An acoustic signal discrimination system of claim 28, wherein the 

2 computing of the output signals is based also on internal states computed on at least 

3 one previous time instant. 

1 35. An acoustic signal discrimination system of claim 28, wherein the 

2 computing of the output signals is/based also on at least one of the input signals, the 

3 output signals, previously received input signals, and previously computed output 

4 signals. 

1 36. An acoustic signal discrimination system according to any one of 

2 claims 28-35, wherein the signal separation architecture is defined by a state space 

3 representation that establishes the relationship between the input signals and the 

4 output signals. 

1 37. An acoustic signal discrimination system according to claim 36, 

2 wherein the computing of the output signals is based also on current states of the state 

3 space architecture. 
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1 38. An acoustic signal discrimination system according to claim 36, 

2 wherein the computing of the output signals is based also on states of the state space 

3 architecture computed on at least one previous time instant. 

1 39. An acoustic signal discrimination system according to claim 36, 

2 wherein the computing of the output signals is based also on the current states and 

3 also the states of the state space architecture computed on at least one previous time 

4 instant 

1 40. An acoustic signal discrimination system according to any one of 

2 claims 36-39, wherein the state space representation is mapped onto a finite impulse 

3 response (FIR) filter. 

1 41 . An acoustic signal discrimination system according to any one of 

2 claims 36-39, wherein the state space representation is mapped onto an infinite 

3 impulse response (IIR) filter. 

1 42. An acoustic signal discrimination system according to any one of 

2 claims 36-39, wherein the state space representation is generalized to a nonlinear time 

3 variant function. 
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